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ABSTRACT 


The  results  of  a  feasibility  study  for  determining 
aerosol  size  distribution  using  the  extinction-sedimentation 
inversion  technique  are  presented.  It  was  found  that  the 
height  of  the  dissemination  chamber,  may  affect  the  value  of 
the  results. 
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1.0  INTRODUCTION 


Particle  size  is  a  factor  which  determines  the  ability 
of  an  aerosol  to  diffuse  light.  An  aerosol's  ability  to 
diffuse  light  will  vary  according  to  its  particle-size 
distribution. 

Several  methods  exist  for  measuring  the  particle-size 
distribution  of  an  aerosol.  First  of  all,  there  are  optical 
methods  using  a  wide  variety  of  instruments  from  the 
microscope  to  the  Doppler  effect  laser  velocimeter.  There 
are  also  non-optical  methods  such  as  the  use  of  fine  grids, 
cascade  impactors,  aerosol  sedimentation  as  a  function  of 
time,  and  the  Model  TA  II,  Coulter  counter  made  by  Coulter 
Electronics  Inc. 

This  report  presents  the  results  of  a  feasibility  study 
for  setting  up  an  optical  method  based  on  measuring  the 
extinction  coefficient  to  determine  the  particle-size 
distribution  of  aerosols  undergoing  sedimentation  in  the 
large  DREV  silo-shaped  chamber. 

This  measurement  technique,  known  as  ESIT 
(Sedimentation-extinction  Inversion  Technique)  is  applied  to 
an  aerosol,  powdered  glass  No.  3419-2  made  by  Ferro 
Industrial  Products,  Ltd.,  to  determine  its  particle-size 
distribution.  In  order  to  check  whether  the  validity 
conditions  of  this  technique  were'  met,  a  study  was  made  on 


the  mass  of  aerosol  deposited  by  sedimentation  on  filter 
paper  as  a  function  of  time. 

The  theoretical  aspects  of  the  ESIT  technique  and  the 
evaluation  of  sediments  on  filter  paper  are  described  in 
Chapters  2  and  3.  Finally,  the  test  results  obtained  are 
presented  in  Chapter  4.  Appendix  A  gives  information  on 
fitting  a  mathematical  function  to  the  measured  optical 
density  values.  Appendix  B  deals  with  computerization  of 
the  study. 

This  work  was  performed  at  the  DREV  between  March  and 
August  1981  as  part  of  NCP  21B05  "Aerosols." 

2.0  OPTICAL  transmission  of  an  aerosol  undergoing 
SEDIMENTATION 

The  theoretical  aspect  of  the  technique  of 
particle-size  distribution  determination  is  described  in 
this  chapter.  The  rate  of  particle  sedimentation  under  the 
influence  of  gravity  is  discussed  first,  followed  by 
transmission  of  a  monodisperse  aerosol  undergoing 
sedimentation,  and  finally,  an  analytical  expression  of 
particle-size  distribution  for  polydisperse  aerosols  is 
developed. 

2 . 1  Sedimentation 

A  spherical  particle  in  an  undisturbed  viscous  medium 
(absence  of  temperature  gradient,  wind,  etc.)  falls  under 


the  influence  of  the  gravitational  field  at  a  constant  rate 
(ref.  1)  equal  to 


v 
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where  r  is  the  radius  of  the  particle  in  cm. 
pis  the  fluid  density  g/cm3. 
fyis  the  particle  density  g/cm3. 

3  is  equal  to  981  cm/sa. 

^tis  the  viscosity  of  the  fluid  g/cm.s. 

with 


(1) 


C 
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1  +  7  (A1  +  A2 


e-A32r/*j 
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where  (/  is  the  mean  free  path  of  the  molecules  in  the  fluid 
Ax  =  1.257, 

A2  =  0.400, 

A3  -  0.55. 

Table  1  gives  the  value  of  coefficient  C  and  the 
sedimentation  rate  of  particles  with  different  diameters. 


TABLE  I 


TRANSPORTATION  PROPERTIES  OF  AEROSOLS  (Ref.  1) 


SPHERICAL  PARTICLES  IN  AIR  AT  20 °C,  1  atm 


a  function  of  their  sedimentation  rate,  namely: 


2 . 2  Optical  Transmission  Through  A  Monodisperse  Aerosol  in 
Particles  in  Process  of  Sedimentation 

A  monodisperse  aerosol  is  an  aerosol  whose  particles 
are  all  of  the  same  diameter  and  hence  fall  at  the  same 
rate.  Figure  1  represents  optical  transmission  through  a 
monodisperse  aerosol  of  uniform  concentration.  At  a  time 
tQ ,  the  aerosol  is  disseminated  instantly.  After  a  time 
fcl  “  tQ'  there  Is  no  longer  any  aerosol  above  a  height 
H-h  (h  and  h  are  defined  in  Figure  3) .  The  sedimentation  is 
given  by: 


.  *» 

8  C1  ~  c0  (4) 

Substituting  the  value  of  vg  in  Equation  2,  we  obtain 
the  radius  of  the  particles  of  which  the  aerosol  is 
composed,  namely: 


un  i  t*'  orb  I  tra  I  r* 


Figure  1.  Transmission  as  a  function  of  time  of  a  monidisperse 
aerosol  undergoing  sedimentation. 

(a)  Transmission  (%) ,  arbitrary  value.  (b)  Arbitrary  unit  time 
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Now  that  we  know  the  ratio  I  .  /I  ,  where  I  . 

min  o  mm 

and  I  are  the  light  intensities  in  the  presence  and  in 
the  absence  of  aerosol,  we  can  use  Beer's  law  to  find  out 
the  number  of  particles  (per  unit  volume)  which  have  been 
disseminated. 


T  ujisH- 

I  -  I  e  oQ  T  *  In  I/I 
o  o 


(6) 


where  (ref.  1) 


l. 

T  =•  n  ff  K(x  ,  n)  r“-r.(r)dr  dl 

0  0  n  w ) 

where  K  (x  ,  n)  is  the  effectiveness  of  diffusion, 
n 

n  is  the  index  of  refraction, 
xn  is  equal  to  2nr/X  f 

A  is  the  wavelength  of  the  incident  light, 
n(r)  is  the  distribution  function, 
r  is  the  radius  of  the  particle, 

L  is  the  length  of  the  optical  path. 

The  function  K (x^,  n)  cannot  be  evaluated  with  Mie's 
theory.  It  is  of  the  oscillating  type  and  tends  toward  2 

I 

for  high  values  of  x^. 


For  a  monodisperse  aerosol,  the  distribution  function 
can  be  written  with  the  aid  of  a  Dirac  function  centered  on 
the  radius  r^  of  the  particles: 


n(r)  -  N6 (r  -  r) 
P 


(8) 


where  N  is  the  number  of  particles  per  unit  volume. 


After  integration,  and  isolating  N,  it  appears  that 


„  T  _  T 

N  “  Lnr*  K(xn,  n)  2 im^L  (9) 

for  high  values  of  xfl. 

Hence,  it  is  possible  to  determine  the  size  and  number 
of  particles  of  a  monodisperse  aerosol  by  studying  its 
sedimentation  rate  and  the  optical  transmission  through  the 
aerosol . 

2 . 3  Optical  Transmission  through  a  Polydisperse  Aerosol  in 
Process  o f  S e dime ntatio n 


An  analysis  of  the  foregoing  type  can  be  done  for  a 

* 

polydisperse  aerosol  by  using  summations.  However,  Deepak 


and  Vaughan  (Ref.  2)  developed  a  more  analytical  method 
which  is  relatively  simple  to  use  and  gives  the  distribution 
function  directly. 


2.3.1  General  Form  of  Distribution  Function 


at  a  time  t  +  At,  we  have: 


r2(t  +  Ac) 

T(t  +  At)  -  nL  /  K(x  ,  n)  r2  n(r)  dr 
rl  n 


(13) 


Using  the  mean  value  theorum  to  evaluate  these  two 
integrals  and  the  definition  of  the  derivative,  we  obtain: 


dT(t) 

dt 


llm  T(t  +  At)  -  T(t) 
At+0  At 


(14) 


-  xLK(x(t),  n)  n(r(t))  r2(t)  (15) 


From  the  latter  equation,  the  distribution  function  can 
be  isolated: 


.  41(0.  /  *L  K(x(t),  n)  r2(t) 

dt  dt 


n(r(t)) 


(16) 
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2.3.2  Equation  of  Distribution  Function 

If  we  are  to  define  the  equation  of  the  distribution 
function  t 16 ]  completely,  we  must  find  out  the  derivatives 
with  respect  to  time  of  the  optical  density  and  critical 
radius  of  the  particles. 

The  choice  of  a  function  to  represent  optical  density 
enables  the  distribution  function  to  be  defined 
analytically.  We  chose  the  following  function  to  fit  the 
experimental  points: 


T(t) 


A(eCt  -  1)  +1) 


(17) 


T0  is  the  value  of  the  optical  density  of  t  =  0,  and 
the  constants  A  and  C  are  determined  with  the  least  squares 
method  after  transformation  of  equation  17  (See  Appendix  A) . 

The  expressions  for  the  derivative  are  obtained  with 
the  aid  of  equations  17  and  10. 


dT(t)  -ACT0eCt 
dt  (A(eCt-l)+l)2 


(18) 


dr(t)  m  l  .-3/2 
dt  2 


The  expression  for  the  distribution  function  becomes: 


n(r(t)) 


2ACT0«Ct  t5/2 


M  Q  \ 


1)+1)2  K(x(t ) ,  n)wLa 


It  is  preferable  to  express  variable  t  as  a  function  of 


r  to  prevent  confusion: 


ACT0e  r2  a 

n(r)  - - - — - - 

(£2 

(A(e  r  -  1)  +  1)2  nLr5  (20) 

For  the  latter  equation,  the  diffusion  effectiveness  K 
has  been  established  at  2.  Equation  20  is  valid  for 
particles  with  radii  greater  than  1  micron. 

3.0  THEORETICAL  MASS  DEPOSITED  AS  A  FUNCTION  OF  TIME 

The  theory  described  in  Chapter  2  which  enables  the 
particle-size  distribution  function  to  be  calculated  assumes 
that  the  particles  undergo  sedimentation  under  ideal 
conditions,  i.e.: 

no  agglomeration  by  diffusion  or  caused  by 
differences  in  sedimentation  rate  of  particles 
with  different  sizes, 
still  air, 

spherical  particles, 

negligible  loss  by  diffusion  on  walls. 

One  way  of  detecting  the  presence  of  processes 
affecting  aerosol  sedimentation  rate,  and  hence  delimiting 
the  range  of  validity  of  the  technique  studied,  consists  of 


measuring  the  actual  mass  deposited  as  a  function  of  time 
and  comparing  it  to  the  calculated  mass. 

The  mass  deposited  by  sedimentation  is  calculated  for  a 
discrete  distribution  and  then  for  continuous  distribution. 

3 . 1  Calculation  of  Mass  Deposited  by  Sedimentation; 

Discrete  Distribution 

Let  n  be  the  number  of  particles  with  radius  r  per  unit 

volume,  let  v  be  the  sedimentation  rate,  t  the  time 
s 

variable,  h  the  height  of  the  aerosol  column  above  the 
filter,  at  time  t  =  0.  and  S  the  surface  area  of  the  filter 
paper. 

If  the  filter  paper  is  empty  at  time  t  =  0,  the  number 
of  particles  deposited  on  the  paper  after  a  time  t  is  given 
by: 

N  -  nSv  t  H(t  -  h/v  ) 

P  8  8  (21) 

The  function  i?  (t  -  h/v  )  is  a  descending  step 
function;  it  has  the  value  1  for  times  shorter  than  h/v 

s 

and  value  0  for  times  greater  than  h/v  .  This  function 

s 

represents  the  fact  that,  after  time  h/v  ,  all  the 

s 

particles  are  deposited. 

If  the  filter  paper  is  empty  at  a  time  t^  the 
expression  for  the  number  of  particles  deposited  after  time 

t  is 


P  “  nSV8  ^  "  tl>  H(t  “  (h/VS  "  tl)) 


If  we  now  consider  particles  with  different  diameters. 


expression  22  becomes 


Np  -  (t  -  M)  nlVl  H(t  -  tt)) 


where  n^  is  the  number  of  particles  per  unit  volume,  and 
is  the  sedimentation  for  particles  of  type  i. 

The  mass  deposited  rather  than  the  number  of  particles 
on  the  filter  paper  is  measured.  The  mass  deposited  is 
equal  to  the  mass  flux  through  surface  A  multiplied  by 
elapsed  time,  namely, 


<t  -  tt)  ^  In  r*  H(t  -  < 


M  1-1  i  1 


-  t  )) 


2ppgr|  1 


where  the  velocity  v^^  has  been  replaced  by  the 


expression 


9  2 

z  PD *r  r  tu  and  r  is  the  particle  radius. 


.2  Calculation  of  Mass  Deposited  by  Sedimentation: 


Continuous  Distribution 


Let  n(r)  be  a  particle-size  distribution  which  is 
continuous  with  respect  to  particle  radius  such  that: 


(25) 


N  -  /  n(r)dr 
0 


where  N  is  the  total  number  of  particles.  The  number  of 
particles  deposited  on  the  filter  paper  at  a  time  t  when  the 
latter  has  been  uncovered  at  a  time  t^  after  the  beginning 
of  sedimentation  is  given  by: 


-  ti)S  /  n(r)  v  (r)  H  (t  -  — ^ 


(r)  "  tl))  dr  (26) 


Hence#  the  mass  deposited  on  the  filter  paper  is: 


8*p  2  - 

V- (t  -  e*> 


7  n(r)r5  H  (t  - 

0  2ppgr 


r-tl))dr  (27) 


Time  t ^  corresponds  to  a  maximum  radius  rM  equal 
cj  l  a 

to  C — at, — ')  f  which  implies  that  at  time  t.  the 

\  p  ~  i  i 

particles  whose  radius  is  greater  than  r^  are  already 

deposited  on  the  plate  covering  the  filter  paper.  At  a  time 

tp,  where  the  filter  paper  is  once  again  "hidden",  all 

particles  with  radius  greater  than  r_  =  — )  ,  but 

r  F 

less  than  rM,  are  deposited.  Their  contribution  to  the 
mass  of  sediment  on  the  filter  is: 


8ir  pp2gS  rM  h 

,(rM>r>rp)"-27“V‘"  [  (7JJ)  "  cl  >  n<r )r&dr 


*  *■« 


*-v- 


i 


where  , 

2  PP8^ 

%<r)  “9m 

Particles  smaller  than  r„  are  deposited  continuously 

£ 

over  the  interval  1 1 1 ,  t?]  .  Their  contribution  to  the 
total  mass  deposited  on  the  filter  paper  is: 


"VO' <‘F  -  I  n(r)rsdr 


The  total  mass  deposited  on  the  filter  is  hence 
written: 


3 . 3  Choice  of  Distribution  Function 

The  mass  deposited  on  the  filter  paper  can  be  evaluated 
if  the  particle-size  distribution  function  is  known.  It 
should  be  noted  that  the  converse  is  also  true.  The 
distribution  function  can  be  determined  by  measuring  the 
mass  deposited  as  a  function  of  time  (ref.  3) . 

Two  distribution  functions  are  presented  in  this 
section.  The  first  was  established  with  a  model  TA  II 
Coulter  counter  (by  Coulter  Electronics,  Inc.)  before 
dissemination  of  the  powder.  The  second  was  measured  with 
the  ESIT  method. 


m 


The  Coulter  counter  measures  the  changing  conductivity 


of  an  electrolyte  placed  between  two  electrodes  when  a 
particle  passes  between  them.  The  signal  recorded  is 
proportional  to  the  volume  of  the  particle.  After  data 
discrimination  and  processing,  the  computer  furnishes  the 
percentages  of  the  mass  of  all  the  particles  analyzed  which 
have  a  diameter  smaller  than  d. 

An  initial  distribution  measurement  (before 
dissemination)  of  the  particles  in  terms  of  the  percentage 
of  the  mass  of  particles  whose  diameter  is  less  than  d  was 
obtained  for  powdered  glass  No.  3419-2.  The  measured 
density  of  these  glass  balls  is  4.2  g/cm3.  The 
particle-size  distribution  in  terms  of  mass  can  be 
represented  by  a  distribution  function  of  the  log 
probability  type  (Figure  2)  . 


(Log  d  -  Log  d  )2 

F(n)  “  m  Jg  0  e  2  Log2 om  < 

m 

where  F (m)  is  the  mass  of  the  particles  with  diameter  d, 
M  is  the  total  mass  of  particles, 
crjjj  is  the  geometric  standard  deviation  by  mass, 
dm  is  the  mean  geometric  diameter  by  mass. 


50%  of  the  particles  have  a  smaller  mass  than  that  of  a 


probability  scale.  (a)  Diameter  (D) ,  microns.  (b) 
Percentage  of  mass  below  D. 


The  curve  in  Figure  2  as  determined  by  the  least 
squares  method  best  represents  the  values  measured  with  the 
Coulter  counter. 


Thus,  using  the  mathematical  expression  representing 
the  best  curve,  or  proceeding  graphically,  we  find  a  mean 
diameter  per  mass  of  10.5  microns.  The  geometric  standard 


TT 
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deviation  per  mass  is  calculated  as  follows  (Ref.  4): 


a  -  _  d  50* 

“  d  50*  d15.87* 


(32) 


where  the  diameters  indicated  correspond  to  15.87%,  50%  and 
84.13%  of  the  mass.  The  mass  is  2.42  for  #3419-2  powdered 
glass. 


3. 3. 2.1  Log  Probability  Type  Distribution  in  Terms  of 
Number  of  Particles 

In  the  preceding  section,  a  log  probability 
distribution  in  terms  of  mass  was  adjusted  to  the 
measurements  of  the  Coulter  counter.  This  distribution  in 
terms  of  mass  must  be  converted  into  distribution  in  terms 
of  number  of  particles. 

The  number  of  particles  with  diameter  d  per  unit  volume 
is  given  by: 


(Log  d  -  Log  d  )2 

P(d)  -  m  Log  .  *  2Los23«  133) 

where  N  is  the  total  number  of  particles  per  unit  volume, 

<r  is  the  geometric  standard  deviation  by  number  of 
particles, 

dg  is  the  mean  geometric  diameter. 
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50%  of  the  particles  have  a  diameter  smaller  than  d^. 

The  geometric  standard  deviation  by  mass  is  the  same 
as  the  standard  deviation  by  number  (Ref.  4)  .  The  mean 
geometric  diameter  by  number  is  linked  to  the  mean 
geometric  by  mass  by  the  relation: 


2 

Log  d  -  Log  d  -  6.9078  Log  o  (34) 

8  ■  8 

Using  the  values  for'  dm  and  cr  of  the  previous 
section,  we  obtain  a  value  of  1.82  microns  for  the  mean 
geometric  diameter  by  number. 

3. 3. 1.3  Relation  between  function  F(d)  and  Function  n(r) 

The  function  F  (d)  cannot  be  used  directly  to  evaluate 
deposits  on  filter  papers.  In  fact,  as  well  as  the  change 
in  variable  from  d  to  r,  we  must  remember  that  a  log 
probability  distribution  is  the  result  of  another  change  in 
variable  (Ref.  5) : 


N _ 

/  2i  Lo  g  0 


(Log  r-Log  r  )2 

Q 

2  Log  a 

8 


(35) 


n(r)dr 


e 


dLog  r 


Evaluation  of  the  mass  deposited  on  the  filter  papers 
calls  for  knowledge  of  the  function  n(r).  The  relative 
proportions  of  the  number  of  particles  was  established  with 
the  Coulter  counter  but  their  density  is  unknown. 

A  mass  of  powder  is  disseminated  in  the  chamber  at 
time  tD  for  a  period  of  30  s.  During  and  after 
dissemination,  during  another  period  of  30  s,  the  inside  of 
the  chamber  is  turbulent  because  the  mixing  fans  are  on. 
Particle  losses  occur  during  mixing  and,  as  a  result,  the 
mass  of  powder  disseminated  is  no  longer  but  rather 

M  .  The  mass  M  is  determined  by  measuring  the  mass 

r  r 

of  aerosol  deposited  by  sedimentation  on  a  filter 
paper  uncovered  at  the  beginning  of  sedimentation. 
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where  S  and  S  are  the  four  areas  of  the  chamber  and  the 
c 

filter  paper  and  the  heights  h  and  H  are  defined  in  Figure 
3. 

The  particle  density  is  given  by 

M  /V 

N  -  p  c  » 

/  G(r  )r3dr  (36) 

J  0 

where  V  is  the  chamber  volume  and 
c 


J 


G<r)  -  1 2n  ilog~Q~ 


(Log  r  -  Log  r  )2 
2  Log2  o 


(37) 


.  Log  e 


3.3.2  Calculation  of  Mass  Deposited  Using  the  Distribution 
Function  Measured  by  the  ESIT  Method 

Calculation  of  the  mass  deposited  on  filter  papers 
using  the  distribution  function  measured  by  the  ESIT  method 
is  relatively  simple.  One  need  only  substitute  the 
distribution  function  n(r)  established  in  equation  20  in 
equations  28  and  29. 

It  should  be  noted  that,  to  calculate  the  mass 


deposited,  it  is  not  necessary  to  determine  the  particle 
density  as  this  is  included  implicitly  in  the  expression  of 
n(r) . 


1 


4.0  MEASUREMENTS  AND  DISCUSSION 


A  particle-size  determination  method  was  presented  in 
Chapter  2  and  in  Chapter  3  evaluation  of  mass  deposits  by 
sedimentation  was  formulated  mathematically.  In  this 
chapter,  the  experimental  and  theoretical  results  are 
presented  and  discussed. 


4.1  Setup  and  Experimental  Method 

Figure  3  illustrates  the  setup  used.  Despite  the  fact 
that  the  ceiling  of  the  silo  is  semi-hemispherical  in  form, 
it  was  considered  that  the  aerosol  was  distributed  uniformly 
in  a  cylinder  11.5  m  high.  An  He-Ne  laser  (632.8  nm,  15  mW 
dc)  by  Spectra  Physics,  and  a  Model  RP313-1  laser  precision 
pyroelectric  detector  whose  aperture  was  reduced  to  0.5  cm 
diameter  were  placed  10  m  from  the  ceiling  of  the  chamber. 
The  detector  was  coupled  to  a  Model  RK  3441,  power 
ratiometer  by  Laser  Precision;  this  instrument  was  connected 
to  a  graphic  recorder. 

Ten  Gelman  filter  papers  with  diameters  4.7  cm  were 
placed  in  cylindrical  boxes  12  cm  deep.  Each  of  the  ten 
boxes  was  covered  with  an  aluminum  plate.  Nine  of  the 
filter  papers  were  uncovered  one  by  one  at  different  times 
t^  with  the  aid  of  strings  attached  to  the  aluminum 
plates.  The  tenth  filter  paper  served  as  a  reference  and  was 
used  to  determine  the  measurement  error  of  the  mass 
deposited. 


V 
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The  aerosol  was  disseminated  with  nitrogen  at  a 
pressure  of  1.03  MPa.  The  aerosol  was  in  a  plastic  bottle 
whose  neck  was  perforated  with  small  holes.  This  was  placed 
10  m  from  the  floor. 

The  sequence  of  events  during  an  experiment  was  the 
following: 

switch  on  mixing  fans  at  t„, 

disseminate  aerosol  for  30  s  at  tQ, 

switch  off  fans  at  time  t  ,  namely  30  s  after 

s 

the  end  of  dissemination 
uncover  filter  papers  at  time  tN, 
evacuate  silo  at  time  t„. 

Since  the  height/diameter  ratio  of  the  boxes  is  high, 
the  deposits  on  the  filter  papers  are  undisturbed  when  the 
silo  is  evacuated  at  time  t„. 

r 

Two  sedimentation  experiments  were  performed  with 
powdered  glass  3419-2  under  the  conditions  described  in 
Section  1.1.  Table  II  presents  the  operating  conditions  for 
the  two  experiments. 
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Figura  3.  Experimental  setup. 

(a)  Expansion  valve,  150  psi.  (b)  Nitrogen  cylinder. 

(o)  Plastic  bottle  perforated  with  holes.  (d)  Mixing 
fan.  (e)  Detector.  (f)  Cover.  (g)  Box.  (h)  Filter 
paper.  (i)  Thumbtack.  (j)  He-Ne  source.  (k)  Graphic 
recorder. 


Table  II 


Mass  dispersed  Mass  collected  Operating  Time  Duration  of 

on  Filter  I  of  mixing  fans  experiment 

after  time 
mixing  fans 
were  switched 
off 


(g> 

(g) 

(s) 

(s) 

El 

1310 

0.0217 

30 

5430 

El 

995 

0.0112 

90 

5430 

4 . 2  Optical  Transmission  and  Distribution  Established  by 
ESIT  Method 

Figure  4  represents  the  optical  transmission  through 
the  powdered  glass  3419-2  undergoing  sedimentation  as  a 
function  of  time  in  the  quantities  given  in  Table  II. 
Figure  5  represents  the  optical  density  in  the  natural  base, 
namely,  -In  I/IQ.  The  curves  drawn  in  continuous  lines 
are  the  best  representing  the  experimental  points  according 
to  mathematical  expression  17.  The  details  of  calculation 
of  the  coefficients  of  the  best  curve  ar  >  given  in  Appendix 
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Figure  5.  Variation  as  a  function  of  time  of  optical  density 
produced  by  dissemination  of  powdered  glass  #3419-2. 

(a)  Optical  density  (-in  I/I  )  (b)  Time  (sec). 

El:  1310  g  disseminated,  mixing  period  after  dissemination: 

.30  sec. 

E2:  995  g  disseminated,  mixing  period  after  dissemination: 

90  sec. 
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Figure  6.  Distribution  functions  obtained  by  ESIT  method  and 
with  Coulter  counter.  _4 

(a)  Distribution  functions  (cm  )  (b)  Particle  diam.  (cm) 

(c)  Coulter  counter  (d)  ESIT 

Els  1310  g  disseminated,  mixing  period  after  dissemination: 
30  sec. 

E2j  995  g  disseminated,  mixing  period  after  dissemination: 

90  sec. 
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Figures  6  and  7  represent,  on  different  scales,  the 
distribution  functions  measured  by  the  ESIT  Method  (Equation 
20)  using  the  mathematical  expressions  of  the  best  curves  of 
optical  density.  The  initial  distribution  functions  of  the 
powders  before  dissemination  determined  by  a  Coulter  counter 
and  calculated  in  Section  3. 3. 1.3  were  plotted  in  continuous 
lines.  It  will  be  noted  that  for  values  with  high 
diameters,  the  distribution  curves  are  parallel  and 
relatively  close  to  each  other.  The  curves  begin  to  have  a 
different  shape  at  10  microns.  It  should  be  noted  that,  for 
diameters  smaller  than  4.6  microns,  the  distribution 
functions  obtained  by  the  ESIT  Method  come  from 
extrapolation  of  the  optical  density  curves. 


Figures  8,  9,  and  10  represent  the  experimental 

sediments  measured  and  those  established  with  the 
distribution  functions  measured  by  the  Coulter  counter  and 
ESIT  Method.  The  curves  of  the  experimental  sediments 
intersect  each  other;  however,  this  intersection  is 
insignificant  because  it  is  within  the  margin  of  error.  The 
curves  of  the  sediments  calculated  with  the  distribution 
functions  established  by  the  ESIT  Method  also  intersect; 
this  is  because  the  distribution  functions  intersect. 

Figures  11  and  12  regroup  the  results  presented  in 
Figures  8,  9, and  10  by  experiment.  The  sediments  obtained 
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Figure  8.  Sediment  on  filters  as  a  function  of  start  of 
exposure,  with  maximum  exposure  time  5430  sec. 

(a)  Time  (sec)  (b)  weight  (g) 

El:  1310  g  disseminated,  mixing  period  after  dissemination: 

30  sec. 

E2:  995  g  disseminated,  mixing  period  after  dissemination: 

90  sec. 


1.0E-01 


l.QE-02 


□  El:  1310  g  dlaa£ain§8,  plrlode  de  atlange  apr£* 
dlaa&ilnatlon:  30  a. 

*  E2:  995  g  diaa£aln6a,  pCrlode  de  ailange  aprta 

dlsa£alnat Ion:  90  a. 


•  0 


1.0E-03 


1 .  DE-OH  1 - ■ — E.A.1  mil - L— 1-1 1 1 H u _ I  1  I  mill  i  I  i  mm - L- 1  i  imp 

X . 0E-01  1 . OE+D0  1.0E*Q1  l.Q£*02  1.0E*03  l.OE'8* 

T«mp«(«  )(c^ 

Figure  9.  Sediment  (calculated  with  Coulter  counter)  on  fil¬ 
ters  as  a  function  of  start  of  exposure,  with  maximum  exposure 
time  5430  sec. 

(a)  Time  (sec) .  (b)  weight  (g) 

Els  1310  g  disseminated,  mixing  period  after  dissemination: 

30  sec. 

E2:  995  g  disseminated,  mixing  period  after  dissemination: 

90  sec. 


in  experiment  El,  and  those  obtained  with  the  aid  of  the 
distribution  function  measured  with  the  Coulter  counter  are 
in  relatively  good  agreement  for  the  filter  papers  uncovered 
at  the  beginning  of  sedimentation.  The  agreement  is  less 
good  for  the  results  of  Experiment  E2.  The  mixing  fans 
operated  for  an  additional  60  s  during  experiment  E2 ,  which 
favored  losses  on  walls  and  particle  agglomeration.  Since 
the  particles  were  larger,  they  fell  more  rapidly  and  hence 
the  quantity  of  sediment  was  larger. 

The  relative  discrepancy  between  the  experimental  and 
theoretical  sediments  (Coulter  counter)  increases  with  time. 
This  is  due  to  agglomeration  by  diffusion  and  agglomeration 
produced  by  the  difference  in  falling  rate  of  particles  with 
different  sizes.  Hence,  the  ratio  between  the  experimental 
and  theoretical  sediments  is  greater  than  2  after  1,000  s. 

Agglomeration  by  diffusion  of  a  monodisperse  aerosol 
can  be  represented  mathematically  by  the  following 
expression: 

<38> 

ut 

where  N  is  the  particle  density,  is  the  collision 
cross-section  and  v  is  the  brownian  motion  of  the  particles. 
The  characteristic  agglomeration  time  by  diffusion  is  given 


Figure  12.  Comparison  of  experimental  and  theoretical  sediments 
(a)  Experiment  E2,  995  g  disseminated,  mixing  period  after 
dissemination  90  sec.  (b)  weight  (g)  (c)  experimental 

sediments  (d)  theoretical  sediments.  Coulter  counter 
(e)  theoretical  sediments,  ESIT  method  (f)  time  (sec) 


Thus,  for  N  =  6.4  x  10 3  cm3,  d  =  1.82  microns,  and 

4 

at  room  temperature,  T  is  approximately  1  x  10  s, 

c 

which  is  a  very  long  period  of  time.  This  leads  to  the 
conclusion  that,  due  to  the  relatively  broad  particle 
distribution  and  the  height  of  the  silo,  agglomeration  by 
diffusion  is  far  less  than  that  caused  by  the  difference  in 
falling  rate  of  particles  with  different  sizes.  With 
reference  to  this  fact,  the  calculations  by  Friedlander, 
(ref.  1)  show  that  the  diffusion  agglomeration  process  is 
dominant  for  particles  with  radii  smaller  than  1  micron  and 
that  agglomeration  caused  by  the  difference  in  falling  rate 
predominates  for  particles  with  diameters  larger  than  1 
micron.  Note  that  these  calculations  are  based  on  the 
interaction  of  particles  with  radii  1  micron  with  particles 
with  radii  varying  between  0.1  and  10  microns. 

4.4  Discussion 

It  was  shown  in  the  previous  section  that  the  measured 
sediments  and  the  calculated  sediments  (Coulter  counter) 
diverge  due  to  particle  agglomeration. 
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Since  agglomeration  is  a  cumulative  process  its  effect 
is  less  significant  at  the  beginning  of  the  experiment. 
This  explains  why  the  results  before  1000  s  agree  relatively 
well  with  the  theoretical  value  of  the  deposits. 

The  ESIT  Method  is  not  valid  if  there  are  interactions 
between  particles.  Hence,  the  results  obtained  with  this 
technique  are  no  longer  valid  after  1000  s.  After  this 
time,  all  the  particles  larger  than  10  microns  are 
deposited.  As  a  result,  the  distribution  functions  obtained 
by  the  ESIT  Method  are  not  valid  for  particles  whose 
diameters  are  smaller  than  10  microns. 

5.0  CONCLUSION 

The  limits  of  the  ESIT  technique  for  determining 
particle  size  were  established  by  comparing  theoretical 
sediments  and  experimental  sediments  on  filter  papers  as  a 
function  of  time.  The  experiments  show  that  the 
considerable  height  of  the  silo  appreciably  increases  the 
particle  sedimentation  time  so  that  agglomeration  becomes 
important.  The  distribution  function  obtained  is  not  valid 
for  particles  smaller  than  10  microns.  It  would  be 
possible  to  determine  the  sizes  of  particles  smaller  than  10 
microns  by  moving  the  laser  and  the  detector  to  the  upper 
end  of  the  dissemination  chamber.  The  height  of  the  aerosol 
column  above  the  measurement  line  would  be  smaller  and  the 
agglomeration  caused  by  the  difference  in  particle  falling 


rate  would  be  less  important.  Moreover ,  since  the 

experiment  would  last  a  shorter  time,  agglomeration  by 
diffusion  would  be  less. 
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TO 


Best  curve  with  shape  Y  = 


VAf 

Af  +  - S. — £ — 

I  Th 

A(e  -1)+1 


We  will  first  attempt  to  linearize  this  equation  and 
then  determine  the  coefficients  giving  the  minimum  standard 
deviation. 


Al.l  Linearization  of  Equation 


A0  "  Af 

Y  -  A  + - ^ - 

f  A(eCt  -  1)  +  1 


[A-l] 


where  t  is  the  time  variable 

Aq  is  the  value  of  y  when  t  =  0, 

A^  is  the  value  of  y  when  t  =  •**, 

A  and  C  are  coefficients  to  be  determined, 


Modifying  [A-l],  we  have: 


v  -  Ao  "  Af  Ct  , 

\  -  y^-  -  A(e  -  1)  +  1 


[A-2] 


Deriving  the  expression  on  the  right 


n 


-  ACe 


Ct 


[A- 3] 


and  taking  the  natural  logarithm  of  the  latter  equation,  we 
obtain: 


Z  =  In  Y+  -  Ct  +  D  od  D  -  In  AC  [A-4] 

which  is  a  linear  equation. 

Al . 2  Calculation  of  Coefficients  for  a  Minimum  Standard 
Deviation 

For  a  linear  equation  of  type  z  =  Ct  +  D,  the 
coefficients  for  a  minimum  standard  deviation  are  calculated 
as  follows  (ref.  6) : 


8qJI  i  -  s^o  82wo  '  8lMl  [A— 5  ] 

2  *  D  "  2 

80  *2  **  ®I  80*2  “  SI 


where 


S0 

S1 

s2 

0 


is  the  number  of  experimental  points, 
is  T.  t±, 
is  £  tj, 
is  2  zL, 
is 


However,  in  the  case  which  concerns  us,  some 


explanations  need  to  be  given. 

The  value  of  coefficient  Aq  is  given  directly  by  the 
coordinates  of  the  first  experimental  point  (0,  ylg) ,  i.e. 
the  optical  density  at  time  tQ .  When  the  time  tends 
toward  infinity,  all  the  particles  are  deposited  and  hence 
coefficient  A^  equals  zero. 

Equations  A-3  and  A-4  suggest  that  first  its  slope  and 
then  its  logarithm  be  taken  to  linearize  the  equation.  In 
concrete  terms,  we  have  equation 


z 


i 


In  A 

0 


i/y 


1+1 


i+i 


[A- 6] 


where  y^  and  y^  +  1  are  experimental  points  taken  at 
t^  and  t^  +  It  should  be  noted  that  since  the  slope 

has  to  be  calculated  the  number  of  z^  is  equal  to  n  -  1 
where  n  is  the  number  of  experimental  points. 
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